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ABSTRACT 

Context. Current theory regarding heavy element nucleosynthesis in metal-poor environments states that the r-process would be 
dominant. The star HD 140283 has been the subject of debate after it appeared in some studies to be dominated by the s-process. 
Aims. We provide an independent measure of the Ba isotope mixture using an extremely high quality spectrum and an extensive x 2 
analysis. 

Methods. We have acquired a very high resolution (R = A/AA = 95 000), very high signal-to-noise (S/N = 1110 around 4554 A, 
as calculated in IRAF) spectrum of HD 140283. We exploit hyperfine splitting of the Ba n 4554 A and 4934 A resonance lines in an 
effort to constrain the isotope ratio in ID LTE. Using the code ATLAS in conjunction with KURUCZ06 model atmospheres we analyse 
93 Fe lines to determine the star's macroturbulence. With this information we construct a grid of Ba synthetic spectra and, using a 
X 2 code, fit these to our observed data to determine the isotopic ratio, f a, which represents the ratio of odd to even isotopes. The 
odd isotopes and 138 Ba are synthesized by the r- and .s-process while the other even isotopes ( 134,136 Ba) are synthesized purely by the 
s-process. We also analyse the Eu lines. 

Results. We set a new upper limit of the rotation of HD 140283 at v sin i < 3.9 km S" 1 , a new upper limit on [Eu/H] < -2.80 and 
abundances [Fe/H] = -2.59 ± 0.09, [Ba/H] = -3.46 ± 0.11. This leads to a new lower limit on [Ba/Eu] > -0.66. We find that, in 
the framework of a ID LTE analysis, the isotopic ratios of Ba in HD 140283 indicate f oM = 0.02 ± 0.06, a purely s-process signature. 
This implies that observations and analysis do not validate currently accepted theory. 

Conclusions. We speculate that a ID code, due to simplifying assumptions, is not adequate when dealing with observations with high 
levels of resolution and signal-to-noise because of the turbulent motions associated with a 3D stellar atmosphere. New approaches to 
analysing isotopic ratios, in particular 3D hydrodynamics, need to be considered when dealing with the levels of detail required to 
properly determine them. However published 3D results exacerbate the disagreement between theory and observation. 

Key words. Stars: individual: HD 140283 - Stars: Population II - Stars: abundances - Galaxy: evolution - Nuclear reactions, nucle- 
osynthesis, abundances 

1. Introduction high, the nuclear reaction 22 Ne(a, n) 25 Mg provides the main 

source of neutrons, however neutro n-capture is mostly weak 
Heavy-element abundances are predominantly due to two s-process (Pignatari & Gallino 2008). The 22 Ne abundance is 
classes of neutron-capture processes, the rapid (r-) process and heavily dependent on the initial CNO abundance and the weak 
the slow (s-) process. For the s-process the beta-decay life- ^-process prod uces little Ba relative to lighter species such as Sr 
| time is shorter than the timescale for neutron-capture. These dGallino et al. 2000). 
two clas ses can be subdived into the main, weak and strong s- The astrophysical origin of the r-process is still relatively 
proce s s (|Clavton & Rassbachll967l : lBusso et alJ200lHThe et alJ unknown. The most widely proposed site for the r-process is 
l2007t ISneden et alJ I2008I) and the main and w e ak r-process when a massive star (M > 8 M^) becomes a c ore-collapse su- 
dTrav aglio et alJ l2004t TWanaio & Ishimarul l2006t llzutani et alJ pernova (IWheeler et al.ll998l : lKaiino et al.l20 02). During a core- 
l2009h . Each neutron-capture process occurs in different en- collapse supernova the neutron flux is believed to be so high that 
vironments. The main s-process occurs in late-type, low- to the neutron-capture timescale is shorter than the beta decay life- 
intermediate-mass stars (1 Mq £ M % 8 Mq), during ther- time. Other possible r- process sites have been co nsidered such as 
mal pulsing on the asymptotic giant branch (AGB). An un- neutron star mergers dFreiburghaus et aO 19991) . however, these 
certain physical event or process is presumed to cause un- seem to have been ruled out as dominant s ources for r-proces s 
processed H to m ix with C-rich m aterial in the He-burning material due to their low rates of occurrence (lArgast et al.l2 004). 
shell to form 13 C dBusso et al ] l2001h . In this environment 13 C Several theoretical scenarios have bee n explored in an effort to 
supplies the necessar y neutrons via the reaction 13 C(a, n) 16 understand this phenomenon ( Wan aio & Ishimarull2006l) . 
(Burbidge et al.lll957l) . In the core He-burning phase of solar- The relative importance of the r- and s-process throughout 
metallicity massive stars, where temperatures are relatively Galactic history depends on the evolutionary timescales of the 
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proposed sites and their elemental composition. The lifetimes 
for massive stars are much shorter than for low- to intermediate- 
mass stars. Typical lifetimes for 25 Mq, 8 Mq, 3 Mq and 
1 M<T) stars are ~ 7 x 10~ 3 , 0.04, 0.4 and lOGyrs respectively 
dRomano et al.ll2005l) . As such, the interstellar medium (ISM) at 
the time at which metal-poor (halo) stars were forming (~ 12 
Gyr ago) should have been enriched by the super novae of mas- 
sive st ars a nd hence , r-pro cess material. Papers bv lSpite & Spitd 
(1978!) and iTruranl d!981l) have been particularly influential in 
est ablishing this f r amew ork, as we now discuss. 

I Spite & Spitd y978) analysed 11 metal-poor halo stars that 
have a metal abundance less than 1/100 the solar metal abun- 
dance. They found that Ba and Y were overly deficient rel- 
ative to Fe. Both elements can be formed via either neutron- 
capture process but are dominated by the s-process in the solar 
system where 81% of Ba and 92% of Y is formed via the s- 
process (Arlandini et al., 1999). The more metal-poor stars in 
Spite & Spite's sample had a greater [Ba/Fefl deficiency than 
[Y/Fe] deficiency meaning that as [Ba/Fe] decreases, [Ba/Y] de- 
creases also. This is because a greater fiuence of neutrons is 
needed to form Ba (Z = 56) than Y (Z = 39) so that [Ba/Y] 
gives a good indica tion of the number of neutrons captured (see 
Seeg er et al] [1965). In contrast, Spite & Spite found that Eu, 
94% of which is formed v ia the r-process in solar system mate- 
rial dArlandini et al.|[l999l) . has the same deficiency as Fe, such 
that [Eu/Fe] remains constant as [Fe/H] increases and is essen- 
tially solar at [Fe/H] > -2.6. 

A consideration of the pos s ible si tes and seed requirements 
for neutron-capture led lTruranl (1 198 lh to postulate that neutron- 
capture-element abundances in metal-poor stars should be dom- 
inated by those synthesized through the r-process. This expec- 
tation arises from the realisation that massive stars are capable 
of producing both the Fe-peak seed nuclei and the high neutron 
fluxes even from very low-metallicity gas, whereas intermediate- 
mass stars, while capable of producing neutrons, cannot pro- 
duce the Fe-peak seed s necessary for the main ^-process (see 
also|Gallino et al. 1998). In Truran's interpretation, the variance 
of [Y/Fe] and [Ba/Fe] with [Fe/H] are exactly what one would 
find if the primary source for nuclei beyond the Fe-peak in 
metal-poor stars is due to r-process nucleosynthesis, while the 
i-process begins to contribute more significantly as [Fe/H] in- 
cre ases giving ri s e to th e increase in [Y/Fe] and [Ba/Fe] seen 
by [Spite & Spite (1978). He reasoned in addition that the en- 
hancement of r-process nuclei would indicate that a prior gener- 
ation of massive stars formed during or before t he formation of 
the Ga laxy. In a more quantitative calculation, Travaglio et alJ 
(1999) examined the metallicity dependence of Ba synthe- 
sized in AGB stars vi a the s-process using the chemical evolu - 
tion model, FRANEC dStraniero et al]ll997tlGallino et al.l[l998l) . 
They found that Ba formed via the ^-process has no significant 
contribution to the Ba abundance in the Galaxy until [Fe/H] 
k, — 1. This supports Truran's hypothesis. 

Different mixtures of odd and even Ba isotopes are produced 
by the r- and s-process. In particular, 134 Ba and 136 Ba are pro- 
duced only by the ^-process due to shielding by 134 Xe and 136 Xe 
in the r-process. Although the spectral lines of different Ba iso- 
topes are not well resolved in stellar spectra, the profile core is 
dominated by the even isotopes while the odd isotopes, which 
experience hyperfine splitting (hfs), have more importance in 
the wings of the line profile, relative to the even isotopes. This 
means that the profiles of Ba n 4554 A and 4934 A are depen- 
dent on the contributions of the two processes. Hyperfine split- 



ting arises from the coupling of nuclear spin with the angular 
momentum of its electrons. The nuclear spin is non-zero in nu- 
clides with odd-N and/or odd-Z, i.e. where the nucleus has an un- 
paired n ucleon. Hyperfine structure has been well documented in 
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isotopes, 

late that the fraction of odd isotope^ of Ba is /odd, j = 0. 1 1 : 
for a pure ^-process mixture of Ba and infer f dd, r = 0.46 : 
for a pure r-process mixture. Their numbers are based on mod- 
els which best reproduce the main s-process using an arithmetic 
average of 1.5 and 3 Mq AGB models at Z = i Zq. The errors 
stated here are our propagation of errors associated with indi- 
vidual isotope abundances for Ba. We show the linear relation- 
ship between r-process contributions (as a percentage) and / 0( jd 
in Fi g. [Ha). 

iMagainl ([[995) attempted to verify Truran's proposal by 
measuring the odd fraction in HD 140283, a well stu died metal- 
poor subgiant at [Fe/H] = -2.5 dAoki et all 120041) . but found 
instead that theory and observations were not comparable. He 
used high-resolution {R = A/AA = 100000) high signal-to-noise 
(S IN ~ 400) data. Magain reported the fractional odd isotope 
ratio, / dd, of Ba to be 0.08±0.06, implying that Ba production 
in HD 140283 is predominantly due to the s-process (see Fig. Q] 
(b)) despite [Fe/H] a nd [Ba/Fe] being very low, [Ba/Fe] = -0.8 
(ISpite & Spitd 1 19781) . The code used to resolve the macro tur- 
bulence and analyse the Ban 4554 A line solves the equations 
of hydrostatic equilibrium under the assumption that the stellar 
atmosphere has a plane-parallel geometry (ID) and local ther- 
modynamic equilibrium (LTE). In contrast some more recent 
analyses, which we describe below, compute hydrodynamical 
3D model atmospheres (3D) where the radiative transfer along 
multiple lines of sight is assessed. 

Th e star was later reanalysed by Lamb ert & Allende Prietol 
(120021) . again assuming ID LTE. They obtained a very high- 
resolution (R = A/AA = 200000) high signal-to-noise (S /N ~ 
550) spectrum about the Ban 4554 A line. They found a value 
for /odd = 0.30 ± 0.21 and concluded that, contrary to Magain's 
result, the star is r-process dominated. A value for /,dd = 
0.30+0.21 would imply an r-process contribution of 54% ±60%. 
We note, however, that the error in their measurement of / dd 
means that their result covers the full range of possibilities from 
a pure s-process mix (/,dd,s = 0.11) to a pure r-process mix 
(/odd.r = 0.46); see Fig. Q] This means that although they state 
that their result indicates that HD 140283 is r-process domi- 
nated, the range of /,dd is too broad to be conclusive. We con- 
sider that an r-process contribution of 54% does not substantially 
imply that the star's neutron-capture elements are dominated by 
those synthesised via the r-process. 

Against this background, we sought to improve the determi- 
nation of the ^-process contribution to Ba in this star to help us 
understand the apparent conflicts. 

While we conducted our study ICollet etalj d2009l) reanal- 
ysed the lLambert & Allende Prieto (2002) spectrum in ID LTE 
and, more significantly, also conducted a new 3D hydrodynami- 
cal analysis of the Ba isotopic fraction of HD 140283. In ID LTE 
they found that / odd = 0.33 ±0.13, meaning that 64% ± 36% 
of the Ba isotopes in HD 140283 are synthesized via the re- 
process. Th e central value (0.33) is l i ttle ch anged from that 
obtained bv lLambert & Allende Prietol (120021) (0.30), which is 
not entirely supris ing since they used the same spectrum, but 
ICollet etafl J2009) quote smaller error bars. This is because 
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Fig. 1. a) Relation between / odd and the r- process contribution calculated fro m lArlandini et alJ d 19991). Coefficients a re given where 
/odd = a x r- process (%) + b.h) LAP02: the lLambert & Allende Prietol d2002l) result for f odd . M95: the[Magain ( 1995) result for / odd . 
CAN ID: the lCollet et ail d2009l) ID LTE result for f ndd . CAN3D: the lCollet et al] d2009h 3D hydrodynamical result for / odd . 



lLambert & Allende Prietol d2002l) adopted the standard deviation 
of macroturbulent b roadening es t imates as the main underlying 
measurement error. ICollet et all {2009) use the standard error 
cr/ yN (where N is the number of Fe lines used) as a measure- 
ment of error. The latter is a more reasonable estimate of the 
error as it is a measure of the uncertainty in the mean estimate of 
the broadening. Their analysis of the line using 3D hydrodynam- 
ics gives a value for / odd = 0.15 + 0.12 meaning only 11% ±34% 
of the isotopes are synthesized via the r-process. This value is in 
good agreement with the solar Ba isotopic mix (/ odd;S5 = 0.16 
implying that only 14% of isotopes formed via the r-process 
dArlandini et al][l999l) ) but is once more at odds with the high 
r-process fraction expected under Truran's hypothesis. 

We have obtained a high resolution (R = A/AA = 95 000) 
very high signal-to-noise (S/N = 870 - 1110) spectrum of 
HD 140283. During the course of this paper we discuss how we 
have constrained the macroturbulence by fitting synthetic spec- 
tra to Fe lines. In a detailed error analysis we show how we 
have improved constraining the macroturbulence, which was a 
major source of error that dominated previous studies that anal- 
yse the Ban 4554 A line in ID LTE. The improvement is partly 
due to the higher quality spectrum we have used in this inves- 
tigation. We also explore the impact of using radial-tangential 
macroturbulence, £ rt, ar, d rotational broadening, v sin i (used by 
ICollet etahl d2009l) ) to help constrain macroscopic broadening. 
We then move on to discuss the method used to re-evaluate the 
r- vs. i-process mix by analysing the Ban 4554 A line and, for 
the first time in this context, the Ban 4934 A line in ID LTE. 
Furthermore we discuss the difficulties in analysing the Ban 
4934 A line due to close blends with other lines. Also because 
of the exceptional quality of the data, we have been able to re- 
vise downward the Eu abundance limit for the star. 

2. Observational data 

Our stellar and ThAr calibration spectra were obtained over 
two nights during the commissioning of the High Dispersion 
Spectrograph (HDS) mounted on the Subaru Telescope. The 



stellar spectrum is the sum of 13 exposures totalling 82 min- 
utes. This gives a S /N = 1110 per 12 mA wide pixel around 
the Ban 4554 A line and a S/N = 870 per 12mA wide pixel 
around the Ba n 4934 A line, as measured from the scatter in the 
continuum of the reduced spectrum. The typical resolution as 
measured from ThAr lines is R = A/AA = 95 000. The spec- 
trum was reduced using a ThAr spectrum to wavelength cali- 
brate the stellar spectrum, with typical RMS errors of 1.5 mA 
dAokietal.1 120041) . We utilise the 4554 A and 4934 A lines as 
both arise from the ground state where hyperfine structure is 
large. Although the 4934 A line is weaker - we measure equiva- 
lent widths W4554 = 20.1 mA and W4934 = 13.6 mA - the hfs of 
the 4934 A line is greater, which means both lines can be useful 
diagnostics. We do not attempt to analyse higher excitation lines 
of Ba which are weaker and have much smaller hfs. 



3. Spectral profiles 

To analyse the two Ban line profiles in our spectrum we 
compared our observed profile to synthetic profiles pro - 
duced by the ID LTE code ATLAS dCottrell & Norrisl 119781) . 
We describe below how Fei and Fen lines are used 
to constrain macroturbulence, and then proceed to anal- 
yse the Ban lines. A ID KURUCZ06 model atmosphere 
(http://kurucz.harvard.edu/grids.html) was used with 
parameters for the star T e g = 5 750 K, [Fe/Hl = -2.5 and mi- 
cro turbulenc^ 2 _£_5l.4Jmi s -1 dAoki et al.ll2004l) and log g = 
3.7 [cgs] dCollet et alJl2009T) . 

3.1. Instrumental Profile 

Two ThAr hollow-cathode-lamp spectra over the intervals 
4102 - 5343 A and 5514 - 6868 A, taken during the observ- 
ing run with the same instrumentation and set-up as the stellar 
exposures used in this study, were used to calculate the instru- 
mental broadening. Using IRAF, the full-width at half-maximum 
(FWHM) and equivalent widths of 993 emission lines were mea- 
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sured. It was found that at a wavelength of 4554 A the ThAr line 
FWHM in velocity space (vi nst ) was 3.31 km s _1 , and at 4934 A 
was 3.25 km s _1 . The error in these measurements is taken as 
the standard error of the mean of the individual measurements, 
cr/ yN, which is 0.01 km s _1 , where cr is the standard devia- 
tion of the individual measurements, which is 0.22 km s _1 . We 
assume here that the ThAr lines are unresolved and hence that 
the measured ThAr line width represents the instrumental broad- 
ening. The instrumental broadening could be slightly less than 
that stated, but the difference is immaterial since, in 33.21 we 
measure the combined instrumental and macroturbulent broad- 
ening without needing to distingu ish between the two contribu- 
tions precisely. lAoki et all (120041) showed that the instrumental 
profile is well approximated by a Gaussian. 

3.2. Macroturbulence 

lLambert & Allende Prietol (120021) established that one of the 
major limiting factors in their analysis was the accuracy with 
which macroturbulence could be measured. They found that 
<5/odd/c>FWHM = -0.51 (kms -1 ) -1 , and hence for ctpwhm = 
0.33 km s _1 they achieved an accuracy in macroturbulence cor- 
responding to <Jf M ~ 0.17, dominating their total error of 0.21. 
It was clear that we would have to improve on this significantly 
to make progress. 

We began to constrain macroturbulence by measuring the 
equivalent widths and the FWHM (in velocity space), v„b s , of 
257 apparently unblended Fe i and Fe n lines by fitting Gaussian 
profiles in IRAF. We used this information to produce Fig. [2] As 
AA/A remains constant with wavelength in an echelle spectrum 
(where AA is the width of the pixel in wavelength units) it is 
possible to use Fig.|2]as a check of the quality of measurements. 
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Fig. 2. FWHM versus equivalent width for the 257 Fe i and Fe n 
measured for HD 140283. 93 lines fall into the range of 10 < W 
(mA) < 50; the mean value for v Q b s in this range is 6.9 km s _1 . 



From Fig. [2] we can see that weaker lines, W < 50 mA, al- 
most form a plateau. Here, v b S remains constant even as the 
Doppler core deepens in lines on the linear part of the curve 
of growth, where the Doppler broadening components are dom- 
inant. At W > 50 mA, pressure broadening become significant 
as the core of the line saturates, so the wings begin to broaden. 
Where W < 10 mA, the uncertainty produced by the finite 



signal-to-noise makes it difficult to measure the lines accurately, 
which is shown by the scatter in this region of Fig. [2] 

Of the 257 Fe lines measured, 93 fell between 10mA< 
W <50mA and were used to constrain macroturbulence (re- 
call that W4554 = 20.1mA and W4934 = 13.6 mA). The aver- 
age value for the observed velocity FWHM, v Q b s , in this range is 
6.9 km s _1 . The full list of measurements can be found in Table 

ED 

3.2.1. Gaussian macroturbulence 

We convolve the synthetic flux spectrum of the star with a 
Gaussian of FWHM v con v which represents the convolution of 
the Gaussian instrumental profile with a Gaussian macroturbu- 
lent profile. For now we assume that the star has no significant 
rotation; we shall return to this point in {33.2.2I Current esti- 
mates of rotation of HP 140283 are vsinz' = 5.0 + 2.0 km s _1 
(Ide Medeiros et al]l2006l) . We create a grid of 385 convolved 
synthetic spectra for 1 1 values of macroturbulence 4.9 km s _1 < 
v conv ^ 6.9 km s _1 in steps Av conv = 0.1 km s _1 and 35 values for 
A(Fefl 4.09 < A(Fe) < 5.45 with steps AA(Fe) = 0.04. Each 
synthetic spectrum covered the wavelength range 4 1 00 - 6900 A 
in intervals of AA = 0.01 A. 

To determine the best fit for v conv we compare our synthetic 
model grid to the observed spectrum employing a y 1 test, y 1 = 
YjiOi-M/) 2 /cr?, where 0, is the observed continuum-normalised 
profile, Mi is the model profile of the line produced using ATLAS 
and cr 2 is the standard deviation of the observed points that 
define the continuum, i.e. cr = (S/N) . A ll 93 Fei and Fen 
lines were individually fitted using ay 2 code dGarcfa Perez et all 
2009). This code allows small wavelength shifts, AA, which we 
discuss below and local renormalisation of the continuum of the 
observed profile for every line. It finds values for AA, A(Fe) and 
macroturbulence that minimize^ 2 for each Fe line analysed over 
a window 0.6 A wide and with continuum windows typically 
0.5 A to 1.0 A on each side of this, depending on neighbouring 
spectral features. Values of v conv found by the^ 2 code for the 93 
lines covering the wavelength range 4118 - 6253 A are shown 
in Fig. [3] The full table of results from the x 2 code is found in 
Table|AJ] 

We use an ordinary least squares (OLS) fit to determine v CO nv 
at the wavelengths 4554 A and 4934 A. The OLS has the equa- 
tion v conv = a/1 + b, where a and b are coefficients of the OLS. 
We find that v conv = 5.75 km s _1 and v conv = 5.76 km s _1 
at the Ba 11 4554 A and 4934 A lines respectively. The error in 
these values, represented by the standard error, is 0.02 km s _1 . 
As the uncertainty in these values is greater than the difference 
between them, we adopted one value for v conv for both Ban 
lines, Vconv = 5.75 ± 0.02 km s . Subtracting the instrumen- 
tal FWHM at 4554 A we find the macroturbulence, vr, to be 
v r = V5.75 2 - 3.31 2 = 4.70 + 0.02 km s" 1 . This value agrees 
well with that found by lAoki et all (|2004). The error in vp is 
given by cr 2 r = cr 2 , + cr 2 , , which is equal to +0.02 km s _1 . 

In using the Fe lines to determine the macroturbulence ap- 
propriate to Ba, it is important that we measure lines forming 
over a similar range of depths in the photosphere. This was 
achieved in the first instance by restricting the equivalent width 
range of the Fe lines to span the two Ba lines (see Fig|2|. In ad- 
dition, we have regressed the v conv measurements against equiv- 
alent width, W, and against excitation energy,^, and find no sta- 
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excitation potential, x, and the equivalent width, W, suggesting 
an astrophysical cause. 

3.2.2. Non-Gaussian symmetric broadening 



A . So far we have adopted a Gaussian macroturbulent broaden- 

7TTr7TTr: ^ :+ «a I ing mechanism. We looked at two other macro-scale broadening 

$ +++ * * + - mechanisms, radial-tangential macroturbulence (f rt ) and rotation 

+ (v sin i) (lGravl2008l Chapter 1 &). Each broadening run was given 

the same atmospheric parameter set; T e ff = 5750 K, log g = 3.7, 

; [Fe/H] = -2.5, £ = 1.4 km s _1 . TableQ]shows the results from 

; fitting the 93 Fe lines using the three broadening types, along 

~ with Gaussian instrumental broadening (Vinst = 3.31 km s _1 ). 

a=(3 .42±3.94)xiq- 5 i The third column indicates how many of the 93 lines were best 

standard error = 0.02 6 = 5.60 + 0.1 9 _ J ~ 

11 fit by that broadening mechanism, as judged by the minimum v 

6000 6500 7000 , „ , ,, 

value tor the three methods. 



5500 
X (A) 



Fig. 3. Values of v conv that satisfy the minimum value for^- 2 for 
the 93 lines {plus symbols). The standard error represents the 
scatter from the mean of each line (cr/ y/N) and a and b are co- 
efficients of the le ast squared fit, y C n nv = aA + b. The Ca i and 
Fe i lines used by lAoki et aTl d2004l) to constrain macroturbu- 
lence have been plotted (triangles) to show the consistency of 
our x 2 fits with their results at red wavelengths. 



tistically significant trend of v conv with W, and a weak (2.5cr) 
trend with x- This suggests that using a stricter restriction on 
the Fe line list would not materially alter the macroturbulent ve- 
locity. In the most extreme case, the value for x = OeV would 
imply v conv = 5.62 + 0.05 km s , which (as we show below) 
would increase / dd by 0.09. 

The x 2 code also determined that A(Fe) = 4.91 +0.01, where 
the error is taken as the standard error. Takin g the solar Fe abun- 
dance to be A(Fe)Q = 7.50 ± 0.05 from iGrevesse & Sauvall 
(119981) . we calculate the metallicity, [Fe/H] = -2.59 ± 0.05, 
where the error in [Fe/H] is the propagation of the statistical 
error in A(Fe)„ and A(Fe)Q but so far excludes the systematic 
errors associated with the imperfect choice of atmospheric pa- 
rameters. That error, based on calculations we provide in $5] 
is around 0.07 dex, giving a total error of 0.09 dex . This is in 
good agreement with metallicity we adopted from lAoki et aTl 
(2004J). We note that there is a n updated list of solar abundances 
given in IGrevesse et al. I d2007h calculated using 3D hydrody- 
n amics, however we decide d to use the ID LTE results given 
in IGrevesse & Sauvall d!998l) as we are working in ID LTE. 

We found that the mean wavelength shift, AA, was -12.0mA 
with a standard deviation cr^i = 3.8 mA. There are several rea- 
sons why we would expect to find a wavelength shift between the 
observed and synthetic profiles. The most likely is an error in the 
approximate radial velocity correction of the star, but line-to-line 
differences require further comment. There could be inaccura- 
cies in the assumed wavelengths in the Fe line list, however the 
Fe line list was produced using the m ost up to d ate data avail- 
able through the IRON PROIECT and iNave etall (1 1994b . where 
wavelengths are quoted to 1 mA, and are believed to be accurate 
to < 1 mA. The RMS error in the wavelength calibration was 
reported as only 1.5 mA dAoki et al J 12004), so the line-to-line 
scatter cr A/) exceeds that error. The excess could be due to the 
inability of ID hydrostatic model atmospheres to model turbu- 
lent motions in a star's hydrodynamic atmosphere. Indeed, the 
residuals shifts were found to depend, at least partially, on the 



Table 1. Comparison of all three broadening types. Column two 
gives the broadening based on all 93 Fe lines for the wavelength 
4554 A, determined using the method described in 33.2.11 The 
errors given are the standard error (cr/ VA 7 ). Column three shows 
how many of the Fe lines were statistically better fits with that 
particular broadening technique. 



Broadening Parameter 


Parameter 


(kms- 1 ) 


# of best fit lines 


v r 


4.70 ± 


0.02 


32 


U 


4.37 + 


0.02 


58 


v sin (' 


3.89 ± 


0.02 


3 



Only three of the 93 Fe lines were fit best by rotational 
broadening, and hence we concluded that using only rotational 
velocities to broaden the lines would be unsound. The derived 
value provides a firm upper limit on rotation, v sin i < 3.9 km s , 
in the case with no macroturbulent broadening. The fact that 
most lines are fit better by a macroturbulent profile empha- 
sises that we have not detected true rotation of the star at this 
3.9 km s _1 level. 

The radial-tangential broadening function within the ID LTE 
framework, provides a better fit than Gaussian macroturbulence 
to almost two thirds of the Fe lines. We present Ba results for 
both macroturbulent prescriptions in 34.41 but unless specified, 
our analysis is conducted using Gaussian fitting. 

4. The Ban resonance lines and the barium isotopic 
ratio 

4.1. Ba ii line structure 

There are five principal, stable Ba isotopes that are formed via 
the two neutron-capture processes. The r- & ^-process produce 
different mixes of odd and even isotopes. The r-process does 
not contribute to two even isotopes, 134 ' 136 Ba, which are pure s- 
process isotopes. The two odd isotopes, 135,137 Ba, and even iso- 
tope 138 Ba are formed from both the s- & r-process. The odd 
isotopes broaden the line and make it asymmetric, whereas the 
even isotopes contribute to the centre of the Ba n line and make 

the core deeper. 

We use the hfs information from IWendt et alJ {l984) and 
IVillemoes et alJ dl993l) to compute energy level splittings for the 
lower and upper levels of 135 137 Ba in the 4554A and 4934A 
lines, and incorporated the isotopic shifts relative to 138 Ba. 
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Table 2. The isotopic and hfs information for both Ba lines. The 
oscillator strengths relative to 138 Ba for each line are given in 
column 3 and the calculated gf-values are given in columns 4 
and 5. 

gf-value 



A ( A) 


Isotope 


Relative strength 


i-process 


r-process 


4553.9980 


137 


0.1562 


0.0210 


0.0471 


4553.9985 


137 


0.1562 


0.0210 


0.0471 


4553.9985 


137 


0.0626 


0.0084 


0.0189 


4554.0010 


135 


0.1562 


0.0049 


0.0594 


4554.0015 


135 


0.1562 


0.0049 


0.0594 


4554.0020 


135 


0.0626 


0.0019 


0.0238 


4554.0316 


134 


1.0000 


0.0429 


0.0000 


4554.0319 


136 


1.0000 


0.1450 


0.0000 


4554.0332 


138 


1.0000 


1.1256 


0.7972 


4554.0474 


135 


0.4376 


0.0136 


0.1663 


4554.0498 


137 


0.4376 


0.0589 


0.1320 


4554.0503 


135 


0.1562 


0.0049 


0.0594 


4554.0513 


135 


0.0311 


0.0010 


0.0118 


4554.0537 


137 


0.1562 


0.0210 


0.0471 


4554.0542 


137 


0.0311 


0.0042 


0.0094 


4934.0288 


137 


0.3125 


0.0197 


0.0441 


4934.0332 


135 


0.3125 


0.0045 


0.0556 


4934.0410 


137 


0.0625 


0.0039 


0.0088 


4934.0439 


135 


0.0625 


0.0009 


0.0111 


4934.0750 


134 


1.0000 


0.0201 


0.0000 


4934.0753 


136 


1.0000 


0.0678 


0.0000 


4934.0768 


138 


1.0000 


0.5265 


0.3729 


4934.0918 


135 


0.3125 


0.0045 


0.0556 


4934.0942 


137 


0.3125 


0.0197 


0.0441 


4934.1025 


135 


0.3125 


0.0045 


0.0556 


4934.1064 


137 


0.3125 


0.0197 


0.0441 



= 4554.0332 A 



o.o 

-40 



-10 
AX (mA) 



Fig. 4. The splitting patterns of the Ban 4554 A line relative to 
138 Ba. The relative strengths of each isotope are normalised to 1 
(Tabled column (3)). 



Importantly, we note that the hfs splitting and isotope shift data 
in these papers is quite simi lar to that of some older studies 
(refer ences can be found in IWendt et al.l [19841 : [Villemoes et alj 
1 19931) . and hence we are confident that there is no significant 
uncertainty in the line structure. The line wavelength structure, 
relative to 138 Ba, is shown in Fig. |H in which each isotope is 
shown with a total strength of 1.0. Using the s- and r-process 



contributions to the five isotopes from Arlandini et al. (1999), 
we construct line lists for pure s-process and pure r-process iso- 
tope mixes for the two Ban lines, adopting total log gf values of 
+0.16 and -0.16 for 4554 A and 4934 A respectively, see Table 
[2] From these two lists we created a further 13 hybrid line lists. 
These covered a range for / dd equal to 0.00 < /o dd < 0-46 where 
A fndd = 0.035. We recognise that according to I Arlandini et al.l 
(1999), the cases with / dd < 0.11 are not achieved astrophysi- 
cally because even a pure s-process mixture has a non-zero con- 
tribution of 135 137 Ba. 



4.2. x test 

The observed continuum was renormalised over a window of 
1 A either side of each of the two Ba lines. A new grid compris- 
ing 23 1 synthetic spectra around each of the two Ban resonance 
lines was produced in ATLAS. Values for v conv and A(Fe), con- 
strained in the last section, were fixed. There were three free pa- 
rameters in the new grid: A(Ba), AA and the r- & s-process con- 
tributions. The^ 2 code allowed small changes in these parame- 
ters exactly like the code described in ^3.2.11 We used 2 1 values 
forA(Ba), -1.40 < A(Ba) < -1.20, where AA(Ba) = 0.01. Each 
synthetic spectrum covered the range 4550 - 4560 A (around 
4554 A), 4930 - 4940 A (around 4934 A) and was computed ev- 
ery 0.01 A. The windows in which both Ba lines were analysed 
was +0.25 A from their centroid. 



4.3. The iron blends at 4934 A. 

It has been documented that the Ban 4934 A line has a known 
blend with a weak Fe i line (ICowlev & Frevl[T989l) . We use the 
inform ation for two Fe lines which are found in iNave et ail 
(1994, their table 2). The relevant data can be found in our Table 
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Fig. 5. Synthetic spectra showing the effect of the Fei blends 
on the 4934 A Ban line. The top plot shows the r-process-only 
isotope fraction and the bottom shows the s-process-only isotope 
fraction. Solid line: the underlying Fei blends. Dashed line: the 
uncontaminated Ban line. Dash-dot line: the overall line profile. 



When these wavelengths are compared with the wavelengths 
for the Ba 4934 A line given in Table|2]it is clear that the two Fe 
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Table 3. Spectroscopic information on the two weak Fe 
lines that are blende d with the Ban 4934 A line as re- 
ported by iNave et ail (119941) . The log gf values reported 
here were supplied by G. Nave (2009, private commu- 
nication) but can be found using the Kurucz database 
http : //www . cf a . harvard . edu/amp/ampdata/kurucz2 3/sekur . html. 



A (A) Ion x (eW) log gf 
4934.0052 Fei 4.15 -0.589 
4934.0839 Fei 3.30 -2.307 



Based on the calculations by lArlandini et al.1 (Il999l) . our 
4554 A result should not be achievable, and corresponds to an 
r-process contribution of -29% (i.e. the s-process contribution is 
equal to 129%). We have also plotted in Fig.|7]the fit and residual 
for the nearest physically possible value for / dd (0.1 1). This fit 
has^ = 7.7. We have also plotted the fit for / D dd = 0.01, which 
is quite similar. 



lines would influence the r-process fraction found by analysis of 
the line. This is shown in Fig.[5]for pure s- and r-process isotope 
ratios. The analysis of the 4934 A line is very sensitive to the 
characteristics of the Fe lines, as we show in 35.21 We include 
these Fe lines in our Ba 4934 A line analysis. 

4.4. The r-process contribution 

The Ba abundances of the two lines are found to be A(Ba)4554 = 
-1.28 and A(Ba)4934 = -1.30. We take the implied Ba abun- 
dance as an average of the two, A(Ba) = -1.29 + . 08. U sing 
the solar abundances calculated in lGrevesse & Sauvall(ll998l) we 
find that for HD 140283, [Ba/H] = -3.46 ± 0.11, and hence 
[Ba/Fe] = -0.87 + 0.14. Errors stated here are calculated in §5\ 
Results from other papers are given in Table [4] It is shown that 
our result for [Ba/Fe] is in good agreement with previous results. 

Our x 1 fitting procedure included a possible wavelength 
shift as a free parameter. The 4554 A line has a wavelength 
shift A/1 = -14.8 mA. The 4934 A line has a wavelength shift 
A/1 = -21.5 mA, possibly because of imperfect modelling of 
the Fe blend in the blue wing. Both lines fall within 3<x of 
the mean wavelength shift found in the Fe lines, -12.0mA with 
cr A/1 = 3.8 mA. 

From the x 2 analysis we find the best statistical fit for the 
4554 A line is / od d = 0.01 ± 0.06. The best statistical fit for 
the 4934 A line indicates a value of / dd = 0.11 + 0.19 (mean- 
ing an r-process contribution of 0% ± 54%). The lcr errors 
stated here arise from uncertainties discussed in $5] and are 
larger for the 4934 A line because of uncertainties associated 
with the underlying Fe blends. Our result is in good agree- 
ment with iMagainl (1 19951) who found that Ba n 4554 A yielded 
/odd _ 0-08 + 0.06, but seems to b e at odds with values found by 
lLambert & Allende Prietol (120021) (/ odd = .30 + 0.21) a nd the 
ID result found for the same spectrum by IColleteT al. (2009) 
(/ odd = 0.33 ±0.13). 

The reduced chi-squared values, x 2 , for the 4554 A line and 
the 4934 A line are 6.6 and 2.0 respectively. The two best statis- 
tical fits for both lines and their residuals (observed - synthetic 
profiles) are shown in Fig. [6] We also plot the synthetic profiles 
for the Ba lines with an r-process contribution of 100%, best 
fit A(Ba) = -1.30 for both lines. It can be seen in the residual 
plots for both lines that the pure r-process fits are very poor. In 
4554 A, x 2 r changes faster with / dd than for the 4934 A line. This 
indicates that although the 4934 A line is broader due to the ef- 
fects of hyperfine structure, the 4554 A line is more sensitive to 
changes in / D dd- This could be both because the 4554 A line is 
stronger (W = 20.1 mA) than the 4934 A line (W = 13.6mA), 
and because the latter has an Fe blend. 




-1.0 E , , , , : 

4553.80 4553.89 4553.98 4554.07 4554.16 4554.25 

A (A) 



Fig. 7. Comparison between the nearest physical fit (solid line) 
where / dd = 0.11 and the best statistical fit (dash-dot line) 
where / dd = 0.01, for the 4554 A line. Also plotted are the split- 
ting patterns for barium relative to barium-138 (see Tabled. 



Table 4. Results from previous studies of HD 140283. 



Paper 


7-eff (K) 


log g (cgs) 


[Fe/H] 


[Ba/H] 


[Ba/Fe] 


SS° 


5727 


3.30 


-2.40 


-3.20 


-0.80 


MMZ* 


5640 


3.10 


-2.73 


-3.86 


-1.13 


GS C 


5690 


3.58 


-2.53 


-3.17 


-0.64 


RNB rf 


5750 


3.40 


-2.54 


-3.45 


-0.91 


MGB" 


5640 


3.65 


-2.30 


-3.10 


-0.80 




5650 


3.40 


-2.40 


-3.43 


-1.03 


MK S 


5650 


3.50 


-2.50 


-3.28 


-0.78 


LAP'' 


5777 


3.74 


-2.70 


-3.79 


-1.09 


CAN' 


5690 


3.67 


-2.50 






GRPA' 


5750 


3.70 


-2.59 


-3.46 


-0.87 



°| Spite& SpkeUl978h. 

b Magainl jl989l) andlZhaq & Magainl dl990h . 
c Gratton & Snedenl dl994T 

rf Rvanetal J (1 19961) . 

e Mashonklna et alj < fl999h . 
f FulbrightH2000h. 

* Mishenina & Kovtvukbl (120011). 

* Lambert & Allende Prietol [ 2002). 
' IColletetal.U2009h . 

' This work. 
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Fig. 6. Panel a) i): The best statistical fit synthetic profile (solid line) for the observed Ban 4554 A line (diamonds) with the residual 
and^ 2 plots below. For comparison, a pure r-process line and residual has been plotted (dash-dot line). The value for A(Ba) has been 
optimised to one that minimises^ 2 , values for A(Fe) and macroturbulence remain the same. Panel a) ii): The^ 2 fit for the 4554 A 
line, the cross shows where the minimum of the fit lies. Also plotted are the splitting patterns for barium relative to barium- 138 (see 
Table|2|. Panels b) i) and ii): Show the same as a) i) and ii) for the 4934 A line. 



When we adopted radial-tangential macroturbulence, it was 
determined that / odd = -0.02 + 0.06 and -0.03 + 0.19 for the 
4554 A and 4934 A lines respectively, with x 2 r = 6.1 and 2.8. 
Errors stated here are assumed to be the same as those calcu- 
lated using a Gaussian macroturbulence as only the broadening 
mechanism differs; the errors in the two broadening techniques 
have the same value. The best fits are shown in Fig. [8] So the 
radial-tangential fit for the 4554 A line is a statistically better 
fit than the Gaussian macroturbulent fit, as seen by the residual 
plots (Figs.[6](a)(i) & (b)(i) vs. Figs. 0(a) & (b)). Both broaden- 
ing mechanisms, which were analysed separately, yield similar 
values for AA, / dd and [Ba/H]. Both indicate a strong ^-process 
signature for barium. Although these / Q dd numbers are again be- 
yond possible physical values, we inform the reader that due to 
the finite confidence in the x 2 r test (discussed further in {0, the 
unphysical values / D dd - 0.01 are not greatly preferred over the 
physical value / D dd = 0.11. 



We have given values for / Q dd for the two Ba lines and we 
have shown that the two lines are in agreement within the stated 
errors. We now discuss those uncertainties and what stellar pa- 
rameters / dd is sensitive to. 



5. Uncertainties and sensitivity tests 

In this section we scrutinise the analysis procedures and statis- 
tical tests employed in this study to determine the likely statis- 
tical and systematic errors. These include errors associated with 
the atmospheric parameters used in constructing the synthetic 
spectra, the calculated macroturbulence and the errors associ- 
ated with the iron lines used in conjunction with the Ba 4934 A 
line. 

Table|5]lists values found for [Ba/H] and / dd by varying the 
temperature and log g of the model atmosphere. There are two 
cases. In case 1 we recalculate [Fe/H] and macroturbulence for 
every value of T e s and log g, whereas in case 2 we fix the macro- 
turbulence and [Fe/H] to values calculated for T e s = 5750 K and 
log g = 3.7. Perhaps the first thing to note from this table is that 
altering temperature by +250 K and log g by +0.3 does not drive 
/odd to an r-process dominated fraction. The errors quoted in this 
paper for HD 140283 are based on uncertainties in T e g and log g 
of +100 K and ±0.1 respectively. 

From Table [5] it is possible to calculate the error associated 
with [Fe/H], A(Ba), and hence [Ba/H] & [Ba/Fe], by examin- 
ing how it is affected by the stellar parameters. We use case 1 to 
calculate these uncertainties. It is shown in Table [5] that gravity 
as very little affect on [Fe/H]. Temperature has a much greater 
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4553.80 4553.89 4553.98 4554.07 4554.16 4554.25 4933. B5 4933.93 4934.01 4934.09 4934.17 4934.25 

A (A) A (A) 

Fig. 8. Panel a): The best statistical fit for the 4554 A line (diamonds) using a radial-tangential velocity profile (solid line). We have 
included a pure r-process, / dd = 0.46, synthetic profile for comparison(afa.?/!-£/of line). Also plotted are the splitting patterns for 
barium relative to barium-138 (see Table|2]i. Panel b): Same as panel a) but for the 4934 A line. 



Table 5. The sensitivity to [Fe/H], [Ba/H], v conv and / dd for different values of T e g and log g. case 1: The sensitivity of all atmo- 
spheric parameters used or calculated in this analysis to changes in temperature and log g. case 2: The sensitivity of [Ba/Fe] and 
/odd to temperature and log g. Macroturbulence and [Fe/H] are fixed to values found for the model atmosphere used in this analysis 
(r eff = 5750, logg = 3.7, Vconv = 5.75 kms" 1 , [Fe/H] = -2.59). Columns (11) and (12) show the sensitivity of the blended 4934 A 
line to changes in the stellar parameters. 



/odd 



r eff (K) 


log g (cgs) 


v conv (kms ') 


[Fe/H] 


[Ba/H] 


4554 A 


4934 A 


(i) 


(2) 


(3) (4) 


(5) (6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 






case 1 case 2 


case 1 case 2 


case 1 


case 2 


case 1 


case 2 


case 1 


case 2 


5500 


3.4 


5.79 


-2.78 


-3.69 


-3.71 


-0.01 


0.02 


0.15 


-0.115 


5500 


3.7 


5.67 


-2.77 


-3.59 


-3.61 


0.05 


-0.01 


0.20 


-0.115 


5500 


4.0 


5.52 


-2.74 


-3.50 


-3.50 


0.11 


-0.07 


0.24 


-0.111 


5750 


3.4 


5.82 


-2.60 


-3.50 


-3.52 


-0.02 


0.04 


0.09 


0.113 


5750 


3.7 


5.75 5.75 


-2.59 -2.59 


-3.42 


-3.42 


0.01 


0.01 


0.11 


0.112 


5750 


4.0 


5.60 


-2.57 


-3.33 


-3.31 


0.09 


-0.03 


0.15 


0.112 


6000 


3.4 


5.84 


-2.43 


-3.34 


-3.34 


-0.03 


0.04 


0.02 


0.252 


6000 


3.7 


5.77 


-2.42 


-3.24 


-3.24 


0.01 


0.02 


0.04 


0.252 


6000 


4.0 


5.69 


-2.40 


-3.15 


-3.14 


0.04 


0.00 


0.05 


0.245 



affect on [Fe/H], altering the ratio by +0.07 dex for every 100K. 
Therefore we find a total uncertainty in [Fe/H] of ±0.07. We 
find that &4(Ba)/(51og g = 0.35. Therefore an error of 0.1 in 
log g implies an error crx(Ba),iog g — 0.04. Similarly we find for 
temperature that <XA(Ba)/#T eff = 0.0007 K -1 . Taking the uncer- 
tainty in temperature to be + 100K we find that cr^^j^ = 0.07. 
Macroturbulence affects the shape of lines but not the equivalent 
width. As such we do not include the uncertainties associated 
with macroturbulence here. Also [Fe/H] has very little affect on 
[Ba/H] when compared to temperature and gravity effects so this 
is not included in our error analysis of [Ba /H], The solar bar- 
ium a bundance is A(Ba)Q = 2.17 + 0.07 dGrevesse & Sauval 
119981) . When these uncertainties are added in quadrature we find 
that [Ba/H] = -3.46 + 0.11. Therefore we find that [Ba/Fe] = 
-0.87 + 0.14. 



5.1. The 4554 A line 

We have stated that potentially the most significant parameter 
that / dd is sensitive to is macroturbulence. The difference be- 
tween columns 9 and 10 in Table [5] show how /,dd is sensi- 
tive to macroturbulence. We find that on average, £>/odd/<^conv - 
-0.7 (km s" 1 )" 1 meaning that for cr Vcmv - 0.02 km s _1 , calcu- 
lated in {33.21 o-f oM sc 0.01. That is, by using a large number of 
Fe lines to constrain v CO nv, we have minimized the impact of this 
error. 

Realistically, when you vary one parameter you alter all 
other parameters to compensate for this change. Increasing tem- 
perature increases the derived macroturbulence, which on its 
own decreases / dd- For case 1 we see that an uncertainty in 
temperature of +100K implies an uncertainty in macroturbu- 
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lence of 0.01 km s to 0.02 km s with increasing gravity (see 
Table [5] column (3)). Using the relation we found between / dd 
and macroturbulence we see that <Tf M = 0.02 to 0.06. In case 
2 (where we only look at how f a && is affected by one stellar 
parameter) we find that for an error in temperature of + 100K, 
cr foM ~ +0.004 to 0.01 depending on log g (Table[5]column (10), 
case 2). 

The uncertainty in log g for HD 140283 is quite small, < 0.1, 
due to its reliable Hipparcos parallax. As gravity affects line 
broadening, we find that log g influences the macroturbulence 
and /odd- Firstly we calculate the effect of log g on macro- 
turbulence (case 1). We find 6v com /6log g 0.4 km s de- 
pending on temperature. So for an uncertainty in log g = 0.1, 
cr Vem s= 0.04 km s _1 . Using the sensitivity we calculated for 
macroturbulence suggests an uncertainty in / Q dd - 0.03. In com- 
parison, the total case 1 sensitivity is (5/ dd/51og g m 0.2 imply- 
ing o-f oM = 0.02. We can see from Table[5]the separate effect that 
log g has on / Q dd when we fix macroturbulence (case 2). Here 
we find that <S/odd/<51og g - -0.17 meaning that an error in log g 
of 0.1 alters / Q dd directly by 0.02. The implication is that some 
of the change in case 1 is driven by the revision of the macro- 
turbulence, and some is driven more directly but in a way that 
partially compensates. 

When examining the effect of microturbulence on / Q dd one 
would expect to see two things. If we allow macroturbulence 
to compensate for the change in microturbulence (case 1) we 
would expect find that / dd is essentially unchanged. If we fix 
macroturbulence and alter microturbulence, / Q dd will change. 
Table [6] shows these two cases. As expected in case 1, the 
macroturbulence is driven up/down when the microturbulence 
is decreased/increased and /,dd is unaffected. In case 2 we see 
the sensitivity in / dd as microturbulence is altered given by 
6f dd/St; = -0.5 (km s -1 ) -1 . Therefore an uncertainty in mi- 
croturbulence of 0.1 km s _1 implies an error in / D dd = 0.05. It is 
case 1 that is relevant to our Ba analysis. 

Table 6. The sensitivity of / Q dd to Temperature and log g are 
fixed at 5750 K and 3.7 respectively, case 1: The sensitivity of 
/odd to microturbulence when macroturbulence is re-evaluated 
to compensate for the change to microturbulence. case 2: The 
sensitivity of /odd to microturbulence when macroturbulence is 
fixed at the value calculated when ^ = 1 .4. 



(kms- 1 ) 


^conv 

(kms" 1 ) 




/odd 




4554 A 


4934 A 




case 1 case 2 


case 1 


case 2 


case 1 case 2 


1.3 


5.81 


0.01 


0.06 


0.11 0.14 


1.4 


5.75 5.75 


0.01 


0.01 


0.11 0.11 


1.5 


5.68 


0.01 


-0.04 


0.11 0.07 



In summary we can assign an uncertainty in / dd 
for the 4554 A line crl = crl + <rl + cr] 

/odd v ™nv r clr log g 

V0.01 2 + 0.02 2 + 0.03 2 = 0.04 (case 1 - remember that v conv 
compensates for any effect ^ has on /odd)- In case 2, where we 
look at the separate effects the stellar parameters have on / D dd, we 
find that for uncertainties in macroturbulence, temperature, log g 
and microturbulence cr /odd = VO.01 2 + 0.01 2 + 0.02 2 + 0.05 2 = 
0.06. Case 1 is probably more applicable, but we adopt the larger 
error, case 2, as a precaution, i.e. +0.06. We now move on to er- 
rors and uncertainties associated with the 4934 A line. 



5.2. The 4934 A line 

In order to assign an uncertainty in / Q dd to the 4934 A line we 
must also explore how uncertainties in the Fe blend (see Table 
[3J affect /odd- 

We explored how the equivalent width of the Fe blend is af- 
fected by temperature and log g. As in Tables [5] and [6] we com- 
puted two cases where we allow macroturbulence and [Fe/H] to 
vary with varying temperature and log g (case 1 ) and where we 
have fixed macroturbulence and [Fe/H] (case 2) - in Table [5] 

The net result of an increase in log g, decrease in macrotur- 
bulence, and increase in [Fe/H], is a small increase in synthe- 
sized Wp e , but these effects are minimal compared to the effects 
that macroturbulence has on / D dd- Consequently we see a similar 
behaviour in / D dd (Table [5] column (1 1)) as that exhibited by the 
4554 A line (column (9)) (a roughly linear increase in / dd with 
log g with A/ od d 4934 comparable to A/ od d,455 4 ). 

In case 2 an increasing temperature decreases the equiva- 
lent widths of the Fe lines. Unlike case 1, macroturbulence and 
A(Fe) are not compensating for the increasing ionisation frac- 
tion meaning that Fe i level populations are decreasing. This de- 
creases the strength of the Fe lines, decreasing their equivalent 
widths. 

As log g is driven up in case 2, we find that the equiva- 
lent widths are decreasing, recall that A(Fe) is fixed in case 2. 
Overall, however, we see little or no change in / dd in column 
(12) in Table 

We also investigate how the log gf values, which are not 
well known for the two Fe lines, affect /odd- The 4934 A line 
is driven to a pure r-process fraction if the Fe blend is elimi- 
nated from the line list, so we would expect that / dd would be 
quite sensitive to log gf. Table [3] shows the parameters of the 
two lines. We analyse the case that the log gf values have an er- 
ror +0.15 as a heuristic estimate. Table|7]shows how /odd for the 
4934 A line is affected by this increase/decrease in Fe strength. 
We have tabulated the results for / dd for all values of tempera- 
ture and log g we use in our sensitivity analysis. It can be seen 
for case 2 at a temperature of 5500 K and where the Fe blend 
strengths have been increased, that / D dd becomes so small and so 
unphysical that our^- 2 program cannot find a minimum solution. 
Similarly this is seen in case 1 at a temperature of 6000 K. It is 
clear from Table [7] that for 4934 A, / dd is more sensitive to the 
uncertainty in the strengths of the Fe lines than the atmospheric 
parameters. We see in case 1 that as we alter the Fe log gf by 
±0.15, / dd is altered by +0.18. This means that when added in 
quadrature to the error discussed in {35.11 we find that for case 1, 
/odd = 0.11 ± 0.18. For case 2, / odd = 0.11 + 0.19. We take the 
error to be the average of the two, so / dd = 0.1 1 ± 0.19. 

5.3. Overall result 

Inverse-variance- weighting the results for 4554 A (0.01) and 
4934 A (0.11) give an overall result / odd = 0.02 + 0.06 when 
macroturbulence is modelled as a Gaussian. When a radial- 
tangential broadening mechanism is used we find that inverse- 
variance-weighting gives an overall result / dd = -0.02 ± 0.06. 

So far the uncertainties discussed in the section have been 
limited to errors in , log g, v conv and log gf. We have not yet 
quantified the impact of finite S /N and possible systematic errors 
associated with a ID LTE analysis. We recall from §13.2. 1 1 that 
systematic errors of order 0.09 may arise from using Fe lines 
to estimate, in ID, the macroturbulent broadening of Ba. We 
discuss this further in $7] We shall now move on and discuss the 
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Table 7. Values of / dd for the 4934 A line found from altering the strength of the Fe lines. By changing the log gf values by +0.15 
we find that / dd alters to compensate for a weakening/strengthening of the Ba wing. We have also included an extra decimal place 
to show that there is a slight change in / dd as we alter the temperature and gravity in case 2. 



T s ff (K) 


log g (cgs) 


Vconv ( km s ') 


[Fe/H] 




[Ba/H 








/odd 


(4934 A) 




(1) 


(2) 


(3) (4) 


(5) (6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 






case 1 case2 


i 

easel case2 


case 


1 


case 2 


case 1 


case 2 






-0.15 


+0.15 


-0.15 


+0.15 


-0.15 


+0.15 


-0.15 


+0.15 


5500 


3.4 


5.79 


-2.78 


-3.67 


-3.73 


-3.71 


-3.81 


0.302 


-0.046 


0.126 




5500 


3.7 


5.67 


-2.77 


-3.58 


-3.63 


-3.60 


-3.69 


0.342 


-0.007 


0.127 




5500 


4.0 


5.52 


-2.74 


-3.48 


-3.53 


-3.49 


-3.57 


0.376 


0.043 


0.126 




5750 


3.4 


5.82 


-2.60 


-3.51 


-3.56 


-3.51 


-3.56 


0.252 


-0.150 


0.275 


-0.134 


5750 


3.7 


5.75 5.75 


-2.59 -2.59 


-3.41 


-3.47 


-3.41 


-3.47 


0.273 


-0.127 


0.273 


-0.127 


5750 


4.0 


5.60 


-2.57 


-3.32 


-3.37 


-3.31 


-3.36 


0.316 


-0.100 


0.265 


-0.115 


6000 


3.4 


5.84 


-2.43 


-3.34 


-3.39 


-3.32 


-3.36 


0.200 




0.373 


0.083 


6000 


3.7 


5.77 


-2.42 


-3.25 


-3.31 


-3.23 


-3.27 


0.221 




0.369 


0.084 


6000 


4.0 


5.69 


-2.40 


-3.15 


-3.21 


-3.13 


-3.17 


0.235 




0.362 


0.083 



europium abundance and the various implications of the Ba and 
Eu results. 



6. Europium abundance limit 

Within our spectral range (41 18-6253 A) there are two Eun res- 
onance lines, 4129.70 A and 4205.05 A. Illonda et al.l d2006l) re- 
port th a t the latte r has a know n blend with a V i lin e. Gifrov etafl 
(119881) . iMagainl (1 19891) and iGratton & Snedenl d 19941) report 
[Eu/H] to be -2.31, -2.49 and -2.41 respectively. However, 
these lines do not appear strongly in our spectrum. This be- 
comes clear when studying Fig. [9] which presents the ob- 
served data and several synthetic spectra. The Eu line lists 
were constructed using hy perfine splitting inf ormation from 
iKrebs & Winkl er ( 1960) and iBecker etail (U993I) . We acknowl- 
edge that more recent hyperf ine splitting information is avail- 



able rrom Lawler et al. (/UUl ), wmcn is in good agreement witn 
IKrebs & Winkled (119601) and Becker et al.l 0993) but we do not 
use that data here. An isotopic ratio of 0.5:0.5 for Eu 151:153 
was chosen for the r- and .^-processes (the solar system is o- 
topic ratio of Eu 151: 153 is 0.48:0.52 (lArlandini et al.|[T9 99)). 
and gf values from iBiemont et al.l d 19821) and iKarner et alj 
(1982) were used. The synthetic spectra were produced using 
KURUCZ06 model atmospheres in conjunction with the ID LTE 
code ATLAS. The macroturbulence and [Fe/H] was set at values 
calculated in $3.2.11 

The Eu 4205 A synthesis includes several blends which 
we have adopted from the Kurucz theoretical database 
(http: //www. cfa. harvard. edu/amp/ampdata/kurucz23/sekur . html). 
The dominant blend is a Vn line at A — 4205.09 A, which 
has a gf value of 0.089 and^- = 2.04 eV. The abundance of V 
was taken as the solar abundance scaled to the metallicity. We 
found no blends associated with the 4129 A line. As these line 
blends are theoretical, we make no claim that the abundances 
we deduce from the Eu analysis are as accurate as the analysis 
conducted on the Ba lines. We interpret the absorption feature 
at 4205.1 A as due to Vn, not Eu, as it is much narrower than 
the synthesised, hfs-broadened Eu line. Moreover if it were 
Eu, not V, it would require an abundance inconsistent with 
the weakness of the Eu 4129 A line. We find that a [Eu/H] 
abundance of -2.80 seems to be a generous upper limit on 
the Eu abundance, rather than a genuine detection and lower 



than the cited detections. Therefore we assign a lower limit 
[Ba/Eu] > -0.66. This marginally excludes an r-process ratio , 
whether we assume an r-process limit set in lBurris et al.l (12 000) 
—0.81 , which were calculate d using the [A nders & Grevessa 
1989) isotopic abundances) or Arland ini et all d!999l) (-0.69). 
A pure ^-process rati o (+1 .45. dBurris et al.1 12000) or +1.13, 
(lArland ini et al. 1 19991) ) or a mixed s- and r-process regime, is 
compatible with the data. Howev er, our [Ba/Eu 1 limit d oes agree 
well with observations found in iFrancois et al.l (120071) for stars 
of similar metallicity to HD 140283. 

We shall now move on and discuss the various implications 
of the results found in this paper and look at possible solutions 
to reduce systematic errors associated with a ID LTE analysis. 



7. Discussion 

We have found that / dd = 0.02 + 0.06, and hence the r-process 
fraction implies a purely s-process signature of Ba in HD 
140283. The [Ba/Eu] ratio, > -0.66, is also marginally incon- 
sistent with a pure r-process regime. The isotope result does not 
entirel y co ntradict prev i ous w ork by lLambert & Allende Prietol 
(120021) and lColletetal.1 d2009l) since, due to the size of their la- 
errors, an s- or r-process isotopic mixture was feasible (see Fig. 
[TJ. Although we find that /odd for the 4554 A line is unphysi- 
cal (at the 1.8<x level) based on lArlandini et al.l {l999), we must 
consider the possibility that the adopted s- and r-process iso- 
tope contributions may not be accurate, as they are based on 
our simplified understandin g of nucleosynthesis, which could 
be flawed. For example, the Arla ndini et al.1 d 19991) calculations 
give a solar-system r-process isotopic ratio, and we cannot be 
certain that this applies in the Galactic halo. However, metal- 
poor stars with r-process enhancements do not least have similar 
neutr on-capture abunda nce patterns to the sun, e.g. CS 22892- 
052 dSneden et al.lfT996h . 

We also question whether the S /N ratio is high enough to 
measure these fractions accurately, and whether a ID LTE anal- 
ysis is an adequate tool in investigating isotopic ratios at these 
high levels of S IN, by looking at the confidence limits of the 
X 1 minima, which we now discuss. The fact that our best fitting 
spectra have^ values significantly greater than 1 (j^ = 6.6 and 
2.0 for the 4554 A and 4934 A line respectively) indicates that 
the x 1 denominator (cr,) is not a good description of the devia- 
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[Eu/H] < - 
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Fig. 9. Synthetic spectra for the Eu 4129 and 4205 A lines for [Eu/H] = -2.80 and for the abundances calculated in three other 
studies of this star (see text). It is clear that they over estimate the strength of both lines. While the 4205 A line includes several 
blends particularly V n 4205.09 A, the 4129 A has none. We show the V n line separately in the right-hand panel (dot-line). 



tion of the model spectrum from the data. We interpret the high 
X 2 r values to indicate that systematic errors are present which 
exceed the random fluctuations in the signal. This influence is 
confirmed by inspection of Figs. [6] [7]&[8] where it can be seen 
that the residuals do not oscillate randomly from one pixel to the 
next but rather seem to meander over a cycle of a few pixels. In 
short, this tells us that the failure of the model profile to match 
to the data exceeds the error due to noise (mostly photon noise) 
in the spectrum, and hence cr,- as judged from the S IN under- 
estimates the true residual. From Fig. [6] it can be seen that the 
best fit under-fits the core of the lines in order to fit the wings of 
the lines better. The r-process contributes more to the wings of 
the line, see Fig. [4] and Table [2] It is interesting to note that we 
under- fit the red wing of the 4554 A line between 4554. 1 1 A and 
455 4.17 A (see F igs. & EJ. Lambert & Allende Prietol (120021) 
and IColletetal.1 (|2009) in ID, see the same residual feature at 
this wavelength in terval. When they reanalysed the line in 3D, 
ICollet etall (2009) appeared to remove the feature in the wing. 
This would suggest that it is a result of convection in a 3D at- 
mosphere rather than a feature induced by inaccuracies when 
calculating the isotopic shifts. We suspect that the error arises 
due to the assumptions used in ID LTE codes that are unable to 
correctly model physical conditions in a 3D atmosphere. 

To explore this further, we have searched for evidence of 
asymmetries in the Fe-line data. We have produced two Fe-line 
plots by co-adding the residuals from all 93 lines to find an av- 
erage residual, shown in Fig. [10] (top panel). The lower panel 
shows the average residual for the 82 Fe lines found to have 
no additional features or close non-Fe lines within the window 
over which the^- 2 analysis is calculated (0.6 A). Lines marked 
with an asterisk in Table lA.fl denote the 1 1 Fe lines that were re- 
jected. For one plot (dash-dot curve) the average residuals (obs- 
syn) are based on synthetic spectra calculated using the aver- 
age wavelength shift (-12.04mA) and average macroturbulence 
(5.75 km s~') for all Fe-lines. The average residual is very asym- 
metric, in that the blue wings are better fit than the red wings. 
The large residuals in the red wings, around 60 - 170 mA from 
the line centre, are clearly not due to errors in the central wave- 
lengths of the Fe lines, which are known to better than 1 mA 




Fig. 10. Average residuals for fits to 93 Fe lines, for two cases. 
It is quite clear that assuming a constant wavelength shift and 
macroturbulence is adequate when working in ID LTE but as- 
suming ID LTE when dealing with high quality data is not. Top 
panel: the reduced-noise residual plot for all 93 lines. Bottom 
panel: the reduced-noise residual plot for 82 lines without other 
close absorption features. Lines not included in lower panel are 
noted with an asterisk in TablelA.il 



We also computed a new set of synthetic spectra using the 
optimised v conv and AA values belonging to each Fe line as found 
from the^ 2 analysis (see Table I A. Q , and recalculated the resid- 
uals (obs-syn) for these. These new residuals were coadded and 
averaged, and plotted in Fig.\lQ\(solid curve). These slightly re- 
duce the overall amplitude of the residuals, as expected since 
we are optimising the fit to each Fe line. However, even after 
allowing each one to be optimised, the residuals are still quite 
asymmetric, with the red wings standing out as having larger 
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residuals than the blue wings. From Fig. [10] we can see that a 
similar red feature to the one we see in Fig.[6]for Ba is present in 
both Fe residual plots, and appears at the same distance from the 
centroid of the Fe lines as the residuals for the Ba 4554 A line 
(~ 100mA). The feature remains whether we optimise v conv and 
AA or not. We suspect this feature may be the result of convec- 
tion in the observed, dynamic atmosphere, similar to the feature 
seen in the Ba spectrum, that underlying assumptions in ID LTE 
do not compensate for 

ICollet et alJ (120091) conduct a 3D analysis of HD 140283. 
Our measurements of the residuals for the Fe lines provide a 
future test of whether 3D modelling produces similar residu- 
als. We hope to explore this at a later date. As the 3D pro- 
cess calculates the velocity field ab initio, there is no concept 
of micro- or macro turbulence in that framework. Consequently, 
ICollet et all §009) ascribe any excess broadening to stellar ro- 
tation, and for this they infer v sin i = 2.5 km s _1 . We note that 
the upper limit on v sin i which we infer for zero macroturbu- 
lence is 3 .9 km s~' . Their va lue is compatible with our limit. We 
note that ICollet etalJ §009) find a lower r-process fraction for 
HD 140283 using a 3D analysis than they find for ID. If our find- 
ings are similar, then this will further accentuate the difference 
between the analysis of HD 140283 and the expectations based 
on Truran's hypothesis. 

One possible alternative explanation of the asymmetries is 
that we are seeing the combined spectra of more than one star, 
offset in velocity. It may be difficult to generate the observed 
levels of asymmetries for a realistic second star, and we have not 
attempted to do so, but note this possibility nonetheless. We also 
note that the radial velocity of HD 140283 has been steady t o 
±0.35 km s _1 over long periods of time dLucatello et alJ l2005). 
decreasing the likelihood that it is a binary. 

Finally we shall move on to discuss the conclusions we 
have drawn from our analysis of the isotopic ratio of Ba in 
HD 140283. 

8. Conclusions 

We have used very high quality data (S IN = 870 - 1110, 
R = A/AA = 95 000) to analyse the Ba isotopic fraction and 
the Eu abundance limit in the metal-poor subgiant HD 140283. 
We obtain [Fe/H] = -2.59 + 0.09, [Ba/Fe] = -0.87 ± 0.14, 
and [Eu/H] < -2.80. Using a ID LTE analysis, we find f 0< a = 
0.02 ± 0.06, corresponding to a Ba isotopic fraction which indi- 
cates a 100% contribution by t he s-process. T his result contra- 
dicts the t heory put forward b y Trura n| ( [ 198 lh . The result pub- 
lished bv lLambert & Allende Prietol d2002l) has error bars which 
are too broad to allow one to state conclusively that HD 140283 
is r-process dominated. We have set a new lower limit to the 
[Ba/Eu] ratio, [Ba/Eu] > -0.66. This lower limit marginally 
rules out a pure r-process ratio in HD 140283, consistent with 
the isotopic fraction for barium. A new high resolution spectrum 
with a greater S /N around the 4129 A line is needed to constrain 
a genuine abundance for Eu. 

We have also carried out a careful examination of the 4934 A 
line, which is more sensitive to the effects of hyperfine-splitting. 
We found that, due to the lack of laboratory gf data surrounding 
the Fe blend affecting the wings of this line, it is less effective as 
a tool to analyse the isotopic mixture than Ba 4554 A. 

By examining the spectral residuals for 93 Fe lines and for 
Ba 4554, 4934 A, we find strong line asymmetries in the red 
wing. These may show the shortcomings of using a ID LTE 
analysis to explore isotope ratios; using a more sophisticated 3D 



analysis may be warranted. We are looking to take this work fur- 
ther in th e future and analys e HD 140283 using a 3D code. We 
note that ICollet et alJ d2009l) find a lower r-process fraction for 
HD 140283 using a 3D analysis than they find for ID. If our find- 
ings are similar, then this will further accentuate the difference 
between the analysis of HD 140283 and the expectations based 
on Truran's hypothesis. 

We set a new limit on the rotation of HD 140283: vsinz < 
3.9 km s -1 . 
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Appendix A: Online Material 

Table A.l. The list of iron lines used to constrain v conv including measured W and v b s . 



Measured data from observed spectrum 






Results from^ code 




Wavelength^ (A) 


(mA) 


W° (mA) 


v bs ( km s ') 


Ion 


F ( km s 


- 1 ) A(Fe)d 


AA (mA) 


xl 


4118.54* 


28.8 


32.2 


6.97 


Fei 


5.61 


4.83 


-8.2 


3.3 


4132.90 


14.6 


14.6 


6.59 


Fe i 


5.85 


4.99 


-9.6 


4.6 


4134.68 


26.9 


28.1 


6.84 


Fe i 


5.42 


4.91 


-8.8 


1.6 


4137.00 


11.1 


11.7 


7.07 


Fe i 


5.58 


4.79 


-6.9 


0.9 


4143.41 


34.2 


33.7 


7.25 


Fei 


5.61 


4.83 


-8.3 


5.4 


4147.67 


22.8 


23.5 


6.61 


Fei 


5.52 


4.92 


-11.4 


1.9 


4153.90 


18.2 


19.2 


6.96 


Fei 


5.57 


4.91 


-8.5 


0.7 


4154.50 


21.0 


22.2 


6.65 


Fei 


5.65 


4.84 


-8.6 


2.1 


4154.81 


16.1 


15.9 


6.65 


Fei 


5.40 


4.90 


-9.1 


0.9 


4156.80* 


22.7 


24.6 


7.18 


Fei 


6.00 


4.98 


-8.9 


5.3 


4157.78 


13.7 


13.2 


6.71 


Fei 


5.86 


4.90 


-10.7 


1.4 


4174.91* 


14.3 


14.7 


6.58 


Fei 


5.58 


4.96 


-11.5 


8.7 


4175.64 


20.3 


19.9 


6.63 


Fei 


5.65 


4.97 


-8.9 


4.4 


4176.57 


11.6 


11.9 


6.90 


Fei 


5.90 


5.06 


-9.6 


2.5 


4178.86 


18.5 




6.86 


Fe ii 


5.82 


4.95 


-4.7 


3.2 


4181.76* 


38.8 


39.7 


7.13 


Fei 


5.80 


4.91 


-9.2 


3.4 


4184.89 


17.1 


17.1 


6.86 


Fei 


5.80 


4.89 


-10.5 


1.8 


4187.04 


46.8 


46.4 


7.22 


Fei 


5.79 


4.84 


-11.1 


4.7 


4187.80* 


49.5 


50.2 


7.46 


Fei 


6.07 


4.87 


-12.1 


11.7 


4191.43 


39.5 


38.8 


7.04 


Fei 


5.93 


4.87 


-7.8 


7.2 


4199.09 


48.7 


49.0 


7.32 


Fei 


5.86 


4.80 


-8.7 


14.1 


4210.34 


31.9 


32.2 


7.40 


Fei 


5.84 


4.89 


-9.8 


2.0 


4217.54 


12.0 


13.9 


7.13 


Fei 


5.46 


4.85 


-9.5 


1.2 


4219.36 


23.5 


23.5 


7.25 


Fei 


6.04 


4.92 


-8.6 


2.4 


4222.21 


29.0 


28.5 


6.90 


Fei 


5.88 


4.88 


-8.9 


6.8 


4225.45 


13.2 


14.4 


6.70 


Fei 


5.52 


4.95 


-7.7 


0.9 


4233.17* 


43.0 


43.2 


7.09 


Fe ii 


5.97 


4.87 


-3.8 


9.5 


4233.60 


43.4 


43.7 


7.10 


Fei 


5.81 


4.87 


-10.0 


7.2 


4238.81 


20.4 


21.6 


7.04 


Fei 


6.00 


4.91 


-9.3 


6.5 


4282.40 


46.9 


45.7 


7.12 


Fei 


5.74 


4.84 


-8.8 


10.0 


4337.05 


36.1 


41.9 


7.01 


Fei 


5.71 


4.90 


-14.3 


3.9 


4352.73* 


24.3 


24.4 


6.77 


Fei 


5.88 


4.91 


-10.2 


22.3 


4416.82 


10.9 


11.8 


6.43 


Fe ii 


5.42 


4.97 


-9.1 


0.4 


4430.61 


13.5 


13.2 


6.73 


Fei 


5.77 


5.00 


-9.7 


1.8 


4442.34 


30.3 


31.2 


7.17 


Fei 


5.88 


4.92 


-9.8 


4.1 


4443.19 


11.7 


9.9 


6.64 


Fei 


5.56 


4.85 


-9.6 


0.3 


4447.72 


24.5 


24.0 


6.90 


Fei 


5.85 


4.94 


-10.9 


4.2 


4461.65 


40.9 


42.4 


6.69 


Fei 


5.39 


5.01 


-13.7 


6.4 


4466.55 


29.4 


30.3 


7.01 


Fei 


5.80 


4.88 


-12.7 


3.0 


4489.74 


11.7 


11.6 


6.76 


Fei 


5.56 


4.97 


-14.9 


1.3 


4494.56 


34.4 


34.6 


6.95 


Fei 


5.73 


4.94 


-13.1 


4.3 


4508.28 


16.2 


17.2 


6.88 


Fe ii 


5.97 


5.01 


-10.2 


4.7 


4515.33 


11.8 


13.5 


6.61 


Fe ii 


5.79 


4.96 


-6.0 


1.2 


4520.22 


11.2 


12.8 


6.63 


Fe ii 


5.68 


4.98 


-10.9 


1.3 


4522.63* 


23.3 


24.7 


7.01 


Fe ii 


6.12 


4.88 


-10.8 


52.4 


4531.15 


22.3 


22.6 


6.70 


Fei 


5.54 


4.90 


-10.6 


0.7 


4555.89* 


15.8 


18.3 


6.76 


Fe ii 


5.81 


4.89 


-12.6 


4.3 


4583.83 


37.4 


38.6 


6.98 


Fe ii 


5.76 


5.13 


-11.3 


2.6 


4602.94 


21.6 


21.9 


6.97 


Fei 


5.75 


4.97 


-11.0 


2.9 


4736.77 


14.0 


14.6 


6.89 


Fei 


5.78 


4.99 


-14.2 


1.7 


4871.33 


37.1 




7.25 


Fei 


5.93 


4.84 


-12.2 


5.5 


4872.14 


27.8 




7.12 


Fei 


5.93 


4.87 


-11.6 


10.9 



a Lines marked with an asterisk denote those excluded from the lower panel of Fig.[l0]due to contamination with other lines found ±0.3 A from 
the Fe line centre. 
b Equivalent widths of Fe i and Fe n lines measured in this work. 

Equivalent widths of Fe i and Fe n lines measured by Hos ford et al.l d2009T) and lLambert & Allende Prietol d2002h . 
d A(X) = log 10 ($g) + 12 
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4890.76 


35.9 




7.09 


Fe i 


5.81 


4.86 


-15.1 


10.3 


4918.99 


38.8 




7.11 


Fei 


5.85 


4.88 


-11.0 


15.3 


4938.81 


12.1 




6.80 


Fei 


5.59 


4.89 


-15.0 


2.3 


4994.13 


13.1 


14.0 


6.50 


Fei 


5.58 


4.88 


-12.4 


0.8 


5001.86 


13.3 




7.01 


Fei 


5.91 


4.87 


-11.1 


2.2 


5006.12 


27.5 




6.98 


Fei 


5.68 


4.84 


-13.1 


1.4 


5012.07 


32.8 


32.0 


6.78 


Fei 


5.58 


4.98 


-17.9 


8.8 


5041.07 


13.2 




6.89 


Fei 


6.15 


5.07 


-12.5 


19.1 


5049.82 


22.6 




6.90 


Fei 


5.80 


4.92 


-15.6 


3.6 


5051.64 


22.8 




6.61 


Fe i 


5.60 


4.97 


-15.6 


5.4 


5068.77 


11.2 




6.76 


Fei 


5.33 


4.85 


-17.4 


0.7 


5083.34 


15.3 


16.0 


6.46 


Fei 


5.56 


4.88 


-14.6 


3.8 


5098.70 


10.2 




7.34 


Fei 


5.52 


5.03 


-15.9 


0.6 


5107.45* 


11.6 


12.0 


6.44 


Fe i 


5.51 


5.01 


-20.3 


5.1 


5110.41 


24.1 


25.0 


6.61 


Fei 


5.54 


5.06 


-18.8 


4.3 


5123.72 


11.0 


12.0 


6.40 


Fei 


5.61 


4.99 


-15.1 


2.7 


5142.93* 


11.8 




7.25 


Fe i 


6.20 


4.95 


-14.5 


46.7 


5162.27 


12.6 




7.25 


Fei 


5.76 


5.06 


-10.3 


0.6 


5166.28 


10.3 




6.33 


Fei 


5.41 


4.97 


-15.4 


0.9 


5171.60 


41.2 




6.86 


Fe i 


5.44 


4.67 


0.3 


18.0 


5191.46 


21.9 




6.86 


Fei 


5.81 


4.86 


-16.5 


10.5 


5192.34 


28.4 




6.85 


Fe i 


5.78 


4.84 


-16.9 


29.3 


5194.94 


24.4 


26.0 


6.67 


Fe i 


5.79 


4.93 


-19.2 


7.2 


5216.27 


19.6 




6.65 


Fei 


5.60 


4.89 


-14.9 


4.7 


5234.64 


12.9 


13.0 


6.57 


Fe ii 


5.60 


4.89 


-15.8 


1.9 


5266.56 


32.9 




7.09 


Fei 


5.87 


4.85 


-14.0 


8.6 


5281.79 


13.5 




6.81 


Fe i 


5.77 


4.86 


-17.4 


5.1 


5283.62 


19.5 




7.02 


Fe i 


5.89 


4.94 


-17.1 


12.1 


5302.31 


10.2 




6.85 


Fei 


5.93 


4.86 


-13.8 


2.8 


5324.19 


35.4 




7.15 


Fei 


5.91 


5.10 


-15.2 


21.5 


5339.93 


12.7 




6.90 


Fei 


5.88 


4.95 


-18.5 


4.7 


5367.47 


11.3 




7.04 


Fei 


5.81 


4.84 


-12.5 


5.0 


5369.96 


14.1 




7.11 


Fe i 


5.75 


4.81 


-13.0 


5.1 


5383.37 


18.4 




7.01 


Fei 


5.65 


4.78 


-15.3 


10.0 


5393.17 


12.2 




7.20 


Fei 


6.08 


4.88 


-14.2 


4.1 


5404.15 


18.4 




7.52 


Fei 


5.65 


4.93 


-4.3 


4.6 


5415.20 


17.5 




7.26 


Fei 


5.80 


5.28 


-13.6 


2.3 


5569.62 


12.9 




7.19 


Fei 


5.97 


4.84 


-15.2 


4.2 


5572.85 


19.7 




7.05 


Fei 


5.94 


4.29 


-16.0 


4.1 


5615.66 


34.1 




7.18 


Fei 


5.96 


4.80 


-15.4 


5.6 


6252.56 


11.3 


14.0 


6.74 


Fei 


5.79 


5.00 


-21.1 


2.6 



